Fast cooling of trapped ions using the dynamical Stark shift gate 
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A laser cooling scheme for trapped ions is presented which is based on the fast dynamical Stark 
shift gate, described in [Jonathan et al, Phys. Rev. A62, 042307]. Since this cooling method does 
not contain an off resonant carrier transition, low final temperatures are achieved even in traveling 
wave light field. The proposed method may operate in either pulsed or continuous mode and is also 
suitable for ion traps using microwave addressing in strong magnetic field gradients. 

PACS numbers: 



Introduction - The ability to laser cool ions to within 
the vicinity of the motional ground state is a key factor 
in the realization of efficient quantum computation 0. 
Various cooling schemes have been suggested and im- 
plemented, achieving lower and lower temperature and 
increased cooling rates. The schemes range from the 
simplest Doppler cooling to more sophisticated 
ion trap cooling methods, including sideband cooling for 
two level atoms 0, 0, 0; Raman side band cooling 
which may be used for ions and for atoms 0,0 an d cool- 
ing schemes based on electromagnetically induced trans- 
parency (EIT) 0| . The existing schemes may be divided 
into pulsed schemes and continuous schemes. 

The two-level and Raman sideband cooling schemes 
employ Rabi frequencies that are typically an order of 
magnitude lower than the trap frequency. This is nec- 
essary to avoid off-resonant excitations which increase 
the average phonons number (n) and thus lead to heat- 
ing. Hence, these off-resonant excitations increase the 
final temperature of the scheme. To overcome this limit, 
a scheme that employs EIT was introduced in and 
demonstrated in |lOj. This scheme uses interference to 
cancel the undesirable off-resonant carrier transition by 
suppression of absorption. This scheme can achieve low 
temperatures as well as high cooling rates, but the exact 
calibration of the laser intensity and the high intensities 
which are needed, makes the scheme challenging. 

In this work we present an alternative approach to this 
problem that relies on the fast Stark shift gate 0, • 
In this scheme off-resonant carrier transitions that would 
lead to heating are forbidden even when operating in the 
traveling wave regime at Rabi frequencies that are of the 
order of the trap frequency. The higher Rabi-frequency 
and the cancelation of the carrier transition promises in- 
creased cooling rates and lower final temperatures. This 
advantage applies both for the pulsed scheme and the 
continuous cooling schemes. It is important to note that 
unlike the original Stark shift gate which requires care- 
ful intensity stabilization to achieve high fidelity [ill 
the Stark shift cooling method is only weakly affected 
by laser intensity fluctuations. Our scheme is applicable 
for trapped ions or atoms and is also suitable for an ion 



trap scheme employing microwave addressing method in 
strong magnetic field gradients as proposed in [lflj . The 
presented scheme satisfies the dark state condition [li). 
i.e. the ions are pumped into a zero-phonon state that 
no longer interact with the the applied laser fields. 

This paper is organized as follows, We first describe 
the operation of the Stark shift gate and then explain 
the idea of using this gate for cooling both in pulsed and 
continuous operation. To analyze our proposed scheme 
we deduce the rate equation in the Lamb-Dicke limit and 
verify the predictions thus obtained with exact numeri- 
cal simulations. In order to discuss the efficiency of the 
method we compare the proposed cooling scheme with 
previously proposed cooling schemes. In order to demon- 
strate the applicability of the scheme to chains of ions we 
present numerical results for the simultaneous cooling of 
three ions. We finish with a discussion and conclusion. 

The Stark-shift gate - In the following we briefly ex- 
plain the basic mechanism underlying the Stark shift gate 
[ill ll^ l . The Hamilton operator of a single ion driven 
by a traveling-wave laser field of Rabi frequency Q in 
the standard interaction picture with respect to the free 
atom and phonons is given by 

H = fl{o- + exp(^[ae- il/t + a)e ivt ] - iSt) + H.c}, 

where H = 1, 5 is the laser atom detuning, v the trap fre- 
quency and r\ is the Lamb-Dicke parameter |l5j . For S = 
and employing rj <C 1, the exponent may be expanded 
to first order in r\ yielding H = Q(a x + rja y (ae~ wt + 
)). In a further interaction picture with respect 
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to Hq = rio~ x , we obtain Hi = ir]Q,'[e 
e i(2Si+v)t a+a^]+H.c. where a+ = |+)<-|. Setting Q = § 
(see also [li[ll[l3) and neglecting off resonant terms 
we obtain 
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Careful numerical investigations demonstrate that H ss 
represents an excellent approximation to the exact dy- 
namics E3 • It is of particular importance to note 
that the Stark-shift gate does not have off-resonant car- 
rier transitions, despite being operated in the traveling 



wave regime, making it a very accurate and fast quantum 
gate. As the carrier transition represents the dominant 
heating term in laser cooling, this cancellation of the car- 
rier transition suggests the use of the Stark shift gate for 
fast laser cooling. The Stark shift gate requires stabi- 
lization of the Rabi frequency at O = ^. However, in 
laser cooling we are only interested in the preparation 
of a target state - the ground state. This suggests that 
the constraints on the stability of the Rabi frequency are 
much less severe in the application of the Stark shift gate 
to laser cooling than for its gate operation. The verifi- 
cation and exploitation of this idea is the purpose of the 
remainder of this work. 

Stark-shift gate cooling - In the following we will de- 
scribe how the Stark shift gate may be used to imple- 
ment both pulsed and continuous laser cooling schemes. 
We begin with the pulsed scheme as this clarifies the ori- 
gin of the speed advantage in both schemes. We pro- 
ceed in close analogy to the pulsed scheme based on 
the regular gate [||. Assume that the system starts in 
state | l)\n). First a resonant laser pulse implements 
| i)\n) — » |— )\n). Then the Stark shift gate creates the 
rotation |— ) \n) — > |+)|n— 1). A resonant pulse then 
implements |+) \n — 1) — ► |f ) \n — 1) hence moving the 
system to a dissipative level. Spontaneous decay predom- 
inantly leads to |f) \n — 1) — > ||) \n — 1) and the process 
as a whole leads to the net loss of one phonon. This 
set of pulses is repeated until the limiting temperature 
is achieved. This scheme yields lower final temperatures 
and is faster under ideal conditions than the usual side- 
band cooling. It is evident that the bottleneck of the 
scheme is the |— ) |n) — > |+) |n — 1) and it is this process 
that is accelerated by the use of the Stark shift gate. The 
above scheme may also be implemented in the microwave 
regime .13] where instead of the cooling laser, a combina- 
tion of a microwave and a steep magnetic field gradient is 
used and the coupling to the dissipative level is achieved 
by resonant lasers. 

The continuous scheme is most efficiently implemented 
in a three level configuration as shown in Fig. ^ The 
Hamiltonian of this system is given by 

H = u>i \gi) (gi\ +uj 2 \gz) (ffa| \e) (e| +va)a 

- fi (\gi) (e\ e -*("i+ A >* + \g 2 ) (e| e -^+A)* + H.ej 

- n c (\ 9l )(g 2 \e^-^+H.c). (2) 

The laser on the meta-stable \g\) <-> \g 2 ) transition is 
used to implement the Stark shift gate by tuning its Rabi 
frequency to the resonance condition Q c = ^. This im- 
plements the rotation |— ) \n) — > |+) \n — 1). The two 
lasers on the |<?i) «-» |e) transition are chosen to have 
equal Rabi frequency and detuning and thus couple the 
|+) = -^(1^) + \g 2 )) state to the dissipative level e) 
with spontaneous decay rates 2Fi and 2r2 to the states 
\gi),\g 2 ) respectively. The two lasers leave the state — ) 
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FIG. 1: The required atomic level structure is given in the left 
upper corner. A cooling Laser fi c manifests the Stark shift 
gate, creating a rotation — ) \n) <-> +) \n — 1). The fii and 
Q.2 Lasers couple only the |+) state to the dissipative level 
and leave the |— ) state decoupled. The cooling process is 
schematically depicted in the lower right part of the picture. 

decoupled from the dissipative level. It is important to 
note that the lasers on the \gi) <->• |e) do not couple the 
internal and the phonon degrees of freedom, i.e, their 
Lamb-Dicke parameter should be made as small as pos- 
sible by adjusting the angle between the laser and the ion 
trap. 

The basic cooling cycle in the continuous scheme works 
as follows. The laser with Rabi frequency il c creates the 
rotation |— ) |n) <-> |+) \n — 1) based on the Stark shift 
gate. The lasers couple the |+) \n — 1} state to the dis- 
sipative level |e) \n — 1), which dissipates back to one of 
the \g%,n — l)'s which are superpositions of |+) \n — 1) 
and |— ) \n — 1). If the decay is to the — ) \n — 1), one 
phonon has been lost and cooling is achieved, whereas if 
the decay is to the |+) \n — 1) state the cycle is repeated. 
Thus on average phonons are dissipated out of the sys- 
tem. Note that this mechanism is independent of the sign 
of the detuning A. 

Analytical Results - Following the heuristic arguments 
presented so far we will now derive analytic expressions 
for the cooling rate and final temperature of our scheme 
based on a Master equation approach. The dynamics and 
the final state can be obtained from the master equation, 

^- p =l[H,p]+Cp, (3) 

where H is the Hamiltonian eq. J2J and C is the Li- 
ouville operator describing the dissipative channels of 
from level |e) to |<7j). In the following we use the tech- 
niques presented in [TtI Hsl| . In this formalism, the 
equation of motion for the phonons is obtained by adi- 
abatic elimination of the internal degrees of freedom of 
the ions. The final temperature and the cooling rate are 
obtained by expanding the Liouville operator and the 
density matrix in eq. © in the Lamb-Dicke parameter 
C = Co+ridW^+Oir) 3 ), p = p a +j ]Pl +r 1 2 p 2 + 0(r) 3 ). 



3 



Inserting this into eq. results in three equations, 
each for a different power of the Lamb-Dicke param- 
eter as in 17). From these three equations the final 
temperature and the cooling rate are calculated. The 
validity of the expansion in the Lamb-Dicke parameter 
at the point fi c = ^ is restricted by the three condi- 
tions Ti/rj, v 2 r), Avr\ <C ft 2 - The result of this approach 
is a rate equation for the phonon probability distribution 
P(n) 

—P(n) = rf(A\{n + l)P{n+l) - nP(n)} 
at V 

+A+[nP(n - 1) - (n + l)P(n)]) (4) 



where 
A-(u) = 



2rn 2 n 2 



t 2 {v - 2n c ) 2 + (2n 2 + 2fi c )(A -v + n c )y 

(5) 

and A + = A_(— From these expressions the final 
population and the rate are deduced to be 



(n) = 



-4. 



A- - Aj 



(6) 



T 2 {v - 2n c ) 2 + {2VL 2 + {v- 2fi c )(A -v + il c )) 2 
~ 4v (2QJ? 2 + (A + 3ft c ) 2 + 2 (n c (A + fij - ft 2 ))) 
TT=r ? 2 (A_-A + ), (7) 

where VF is the cooling rate. Numerics shown in Fig. |3| 
confirms that the (n) is minimized at the Stark shift gate 
point fl c « f . At Q c = | we find IF = irr/Vfi 2 (fi- 4 - 
4[4r 2 z/ 2 + (3i/ 2 + 2Aj/ - 2ft 2 ) 2 ]- 1 ). The choice ^r/ = 1 
realizes a cooling scheme with a rate W = |?7ft c which 
is close to the inverse of the gate time T = Note 
that this choice saturates the conditions of validity in the 
above derivation, If one of the other conditions is satu- 
rated as well the final population is smaller by the factor 
of 2, which means that the choice A = T is optimal. Nev- 
ertheless, numerical results shown in Fig. |3 corroborate 
that the cooling time is of the order of the gate time. 
Fig. 13 shows numerically and analytically that for the 
final (n) as a function of ft c we obtain a parabola with a 
minimum around the Stark shift resonance. At this point 
we find 



(n> = 



ft 4 



^r 2 + (A + f) (^ 2 + 2(I(A + |) 



,-ft 2 ))). 

(8) 

Inserting the values of the optimum point, 
Yvj],v 2 r\, Ai/rj = ft 2 , yields a temperature propor- 
tional to rj 2 . For ft = v/2 heating A + is weaker than 

cooling A_ whenever ft 2 < K 41 " + ^ 2 ^^^ 12Au ^ . The 
cooling efficiency is not strongly effected by the fluctu- 
ations of the laser. For example, at the optimal point, 
at r\ = 0.05, for intensity fluctuations of 10 percent, the 
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FIG. 2: The Rate at the optimal point. The squares are yrj, 
when Tc is the time that takes to reduce the population from 
(n) = 1 to (n) = 0.01. The blue line corresponds to half a 
period of a Rabi frequency and the green line to the cooling 
rate calculated from the rate equations W = ^r/Q c 
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FIG. 3: (n) as a function of Q c . The minimum is at the 
Stark shift resonance. The continuous line is the rate equation 
result. The squares are calculated from a numerical solution 
of the Master equation for rj = 1/20 and the triangles are for 
r\ = 1/10. In the limit 77 — > The numerics coincide with the 
rate equations. The parameters are T = 10, = 1/10,^ = 1 
and A = 



rate and the final temperature change by less then 50 
percent. 

Numerical results - The results deduced from the rate 
equations are correct in the limit 77 — > 0. In order to 
check the results for finite r\ and see the dynamics, we 
compare these results to a numerical solution of the Mas- 
ter equation. In Fig. 0a plot of the analytical final mean 
phonon number is compared to numerical results and as 
7/ — > the rate equation results are reproduced. In order 
to compare the rates found for the final stages of cooling 
with the actual cooling time we reproduce the dynam- 
ics using Monte - Carlo simulation for the cooling of one 
phonon. Fig. 0] shows the dynamics of the mean phonon 
number. 

The mean phonon number in sideband cooling is (n) = 

(a + 41 (j:) , where a is a geometric factor of the order 
1 |2oL l2l|. In order for sideband cooling to operate at 
the same rate as the proposed cooling, T has to be of the 



4 




mode:1 ,n:1 



mode:1 ,n:2 



27T/V t 



300 



FIG. 4: (n) as a function of time using Monte Carlo simu- 
lation. The parameters used for the simulation: A = 0,r = 
Q,v = l,n c = l/2,fii = fi 2 = 1,?? = 1/10. For these pa- 
rameters the gate time is -4r — 62.5. It could be seen from 
the numerics that the cooling time is in the same order of 
magnitude as the gate time. 





FIG. 5: Monte-Carlo simulation for the cooling of three ions. 
The Rabi frequency is set to f2 c = + 0.15. For three ions 
the frequencies are i>2 = 



order of vj2r\ which means that the sidebands are not 
resolved, hence, sideband cooling could not be applied. 
The cooling rates and temperatures found in our scheme 
could be achieved by the EIT cooling method [2( but in 
this method the laser power should be at least an order 
of magnitude stronger then in the cooling described in 
this work. 

Multi-mode cooling - Due to the cancelation of the car- 
rier transition it suggests itself that our cooling method 
based on the Stark shift gate could be used to cool several 
modes simultaneously. Fig03 shows Monte-Carlo simula- 
tion for the cooling of three modes in a three ion 
chain. Here we chose the laser Rabi frequency near the 
second mode frequency. The cooling rate achieved in this 
way of course slower than for a single ion, but still faster 
than side band cooling. In a different scheme individual 
ions could be addressed with lasers of different Rabi fre- 
quency each at the optimal working point for a different 
mode. Thus the whole chain could be cooled at times of 
the order the Stark shift gate time. 

Conclusions - We have introduced a method for the 
fast cooling of an ion trap which operates with rates 
of the fast Stark shift gate 01 • The final temperatures 
achieved are proportional to n 2 (recoil energy). The cool- 
ing method is efficient both in continuous and pulsed 
operation, achieving lower temperatures and faster rates 
then ordinary sideband cooling. The gate is also suitable 
for ion traps with microwave frequency addressing |l.'i| . 
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